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ABSTRACT 


This document presents the results of the numerical tests performed in evalu- 
ating the Rational Semiconductor Corporation IMP-16 Orbit Determination 
System. Included herein are descriptions of the tests performed and tabulations 
of the numerical results. This document has been prepared in partial fulfill- 
ment of the requirements of Task 971 of Contract NAS 5-24300. 
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SECTION 1 - INTRODUCTION 


The IMP- 16 Orbit Determination System (ODS) is a software product packaged 
on a network of two National Semiconductor Corporation IMP- 16 microproces- 
sors. The system is capable of performing orbit determination from satellite- 
to-satellite tracldng (SST) data in Applications Technology Satellite (ATS) range 
and range-rate (ATSR) format. The estimation scheme used is a Kalman filter, 
a sequential (recursive) estimator. 

This document evaluates this IMP-16 software/hardware system. Section 2 
provides the numerical results of various component and system tests performed 
on the application software. Section 3 discusses overall system timing and per- 
formance tests. 

The system was developed under Task 885. The validation of the filter was 
carried out during the present Task 971. Although the specific configuration of 
a Tracking and Data Relay Satellite (TDRS) which is tracking the Solar Maximum 
Mission (SMM) is used in the CDS as an example, the system has been developed 
primarily as a demonstration system. This is in accordance with the general 
aim of the task, which is to demonstrate the feasibility of using microproc- 
essors in orbit determination work. 

Components tested included the orbit propagator, the state transition matrix, 
the covariance propagator, and the observation model. Each component test 
was either a comparison with a base run or a numerical approximation to an 
analytical expression (e. g. , difference quotient approximation to a partial 
derivative). Brief descriptions of the mathematics involved in testing each 
component are provided herein. Details concerning mathematical specifications 
may be found in Section 3 of Reference 1. 

The application program system tests used the following as observations; 

IMP-16 simulated data (internal test), Goddard Trajectory Determination Sys- 
tem (GTDS) simulated data, and real data run through the IMP-16 preprocessor. 
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SECTION 2 - NUMERICAL EVALUATION OF THE FILTER 


This section presents the numerical evaluation of the application program. 
Section 2. 1 discusses the Individual component tests; Section 2.2 discusses 
the system tests, 

2.1 COMPONENT TESTS 

The Individual filter components tested were the orbit propagator, the state 
transition matrix, Uie covariance propagator, and the observation model. 

2.1.1 Orbit Propagator 

The accuracy of the Rungc-Kutta orbit propagator within the IMP- 16 ODS was 
checked soon after that component was built. The results were presented In 
a memorandum (Reference 2) that forms part of tlie Task 885 file. The basic 
result was that after 1 revolution of the Solar Maximum Mission (SMM), the 
root-sum-square (rss) error in position was 91 meters, although the rss 
error rose as high as 116 meters during the run. The comparison here was 
made between the IMP results and a run with the GTDS using a 15-by-15 Earth 
field with Sun and Moon and an integration step of 10 seconds. 

2.1.2 State Transition Matrix 

The state transition matrix is used in the filter for propagating the covariance 
matrix. A second-order (in At) Taylor series expansion is used to compute 
the state transition matrix. (See Section 3.2.4 of Reference 1 for a complete 
mathematical description. ) 

Each entry of the state transition matrix is a partial derivative of the form 


A numerlcctl approximation to the (I, j) entry can be made by the difference 
quotient 


AX 




AX 


j,nom 


where AX, ^ ^ ® X, .(t,, + dt) - X, (L) 

I, pert, j l,pert,fO ' l,nom' 0' 

AX, X, (t,, + At) - X, (t,,) 

j,nom j,nom'i) j,nom' 0 

Xi .(t + At) w X^-component at time to + At after the Xj- 

’ ^ component was Initially perturbed 


X 


k,nom 


(t) = Xj;-component at time t , with no initial per- 
turbatlons 


The matrix of numerical difference quotient approximations (Figure 2-1) was 
compared to Uie state transition matrix computed by the program (Figure 2-2) 
for a particular epoch and step size (10 seconds). The matrix In Figure 2-3 
contains the error of the numerical approximations relative to the computed 
IMP state transition matrix entries. That Is, the (I, j) entry In this matrix Is 


<p(l,j) 

where <P(l» j) = (It j) entry in the computed state transition matrix 
NUM(l, j) = numerical approximation to <p(i, j) 

2.1.3 Covariance Propagator 

The state error covariance matrix is propagated between filter updates. In the 
absence of state noise the propagation equation, which can be derived from the 
dynamics model equation, Is as follows: 


P -<p^cp 
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Figure 2-X. Approximations to State Transition Matrix Partials 


1.0000019 
7.6240251 E-5 
3,9818046 E*5 
3.8687535 E-7 
1.5054310 E.6 
7.9669221 E.6 


7,5240261 E*6 
1.0000320 
4.8589670 £.5 
1.5054310 E.5 
6.4199011 E-6 
9.7219767 E-6 


3.0818046 E-6 
4,8689670 E.6 
9.9906596 E-1 
7.9669221 6-6 
9.7219767 E-6 
-6.8057764 E-6 


9.9968416 0,0 0.0 

0.0 9.9968416 0.0 


O.Q 

1,000019 
7.6240261 E.B 
3,9818046 5.6 


O.Q 

7.624025 E.S 
1.0000320 
4,8689670 E.B 


9.9958416 
3,9818046 E 
4,8584670 E 
0.99995596 
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Figure 2-2. State Transition Matrix Computed by Program 
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0.0 

8.400 E-3 

-3.900 E-B 

2,883 E-3 

6,326 E.2 

2.119 E-2 
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0.0 
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0.0 

O.O 

—5,844 E"4 
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4.670 E.3 

2.800 E.5 
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Figure 2-3. Relative Error of State Transition Matrix Approxi- 
mations 
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where P « newly propagated covariance 

(p » state transition matrix 

P * previously updated, propagated, or a t riori covariance matrix 

To verify that the IMP ODS was propagating the covariance correctly, the fol- 
lowing test was performed, A diagonal matrix was entered as the a priori 
covariance matrix (Figure 2-4), This matrix was then propagated for 1 minute 
in six steps of 10 seconds each. The propagation was carried out with both 
the IMP-16 ODS and the Onboard Navigation Paclmge (ONPAC) simulator, a 
research tool used for premission and real-time studies of onboard orbit de- 
termination (Reference 3), The results are shown in Tables 2-1 through 2-3. 
The differences shown in Table 2-3 are relative differences, computed as 
follows; 


Diff/i n IMPliiJIzONaji 
- ON(l,j) 

where Dlff(l,j) - (I, j) relative error 

IMP(l,j) = (I, j) element of the IMP propagated matrix 
ON(l, j) = (I, j) element of the ONPAC propagated matrix 

The Irrgest relative error appearing in Table 2-3 Is -2. 62 percent for the 
2 

term. Most of the errors were considerably smaller. 

2,1.4 Observation Model 

Tests were made to verify the accuracy of the observation model. The SST 
range and range-rate measurements modeled by the IMP-16 ODS were com- 
pared to measurements simulated by the ‘GTDS. The first four comparisons 
were reported In a task memorandum (Reference 2), Those results are re- 
printed here in Table 2-4. 
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Figure 2-4. Diagonal A Priori Covariance Matrix 


Table 2-1, ONPAC Covariance Propagation Results 
(1 Minute) 


X 

(km^) 

Y 

(km^l 

Z 

(km^) 

(km^/seo^) 


■HHII 

0.11491 0 

-0.26f-17 E-3 

0.70631 E*4 

0.14977 E-2 

-0.70100 E-5 

0.22941 E-5 


O.l’aWO 

-0.13270 E-3 

-0.701 10 E-5 

0.15068 E-2 

-0.42038 E-5 



0.67413 E-1 

0.14741 E.B 

-0.27092 E-6 

0.74778 E-3 




0.24961 E-4 

-0.12447 6-6 

0.33858 E-7 





0,25119 E-4 

-0,61434 E-7 






0.12461 E-4 
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Table 2-2. IMP ODS Covariance Propagation Results (1 Minute) 


X 

Ikm^) 

Y 

ikm^) 

Z 

(km^) 


(km^/wc*) 

■HHi 

0.11489 0 

-0.2043B E.3 

0.6B777 6*4 

0.14976 E*2 

-0.69509 e-5 

0,22690 E-5 


0.11526 0 

-0.13034 E*3 

-0,69518 E.5 

0.15007 E-2 

-0.41770 E-5 



0.57407 E-1 

0.14442 E-5 

-0.26781 E-5 

0.74772 E-3 




0.24960 E'4 

-0.12397 e.6 

0.33414 E-7 





0.26120 E-4 

-0.61 104 E-7 






0.12461 E-4 


Table 2-3. Errors of IMP Covariance Propagation Relative to 
ONPAC Propagation 



Table 2-4. Observation Modeling 


TIME 

IHMMSS) 

{OCTOBER 15, 1979) 

RANGE (KILOMETERSI/RANGE RATE (CYCLES PER SECOND) 

ORIGINAL 

IMP 

MODEL 

IMP 

ERROR 

FORTRAN 

MODEL 

FORTRAN 

ERROR 

12000 

72181.1323/ 

72180.121 

1.011 

72180.8463 

0.2860 


6698.80435 

— 


- 

- 

12010 

72173.2425/ 

72172.319/ 

0,923/ 

72172.9876/ 

0.2549/ 


5973.84414 

5988.4654 

5.379 

5974.07284 

0.22870 

12020 

72166.1877/ 

72165.271/ 

0.916/ 

72165,9643/ 

0.2234/ 


5347.96519 

5355.4577 

7.492 

5347.69404 

0.27115 

12030 

72159.9687/ 

72159.065/ 

0,903/ 

72159.7761/ 

0,1926/ 


4721,27864 

4726.5946 

5,316 

4721.99304 



0.71440 


NOTES; 1. THE FIRST RANGE RATE CANNOT BE MODELED; A PREVIOUS MEASUREMENT IS NEEDED 
(NONDESTRUCTMODE). 

2. THE ERRORS ARE ABSOLUTE values. 
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This testing was carried further. The graph of the range residuals (GTDS 
simulated range minus the range modeled by the IMP ODS) is shown in Fig- 
ure 2-5. The residuals were computed for each observation (one every 10 sec- 
onds) for the first 3 minutes and for sampled times thereafter. The resulting 
error curve shows that there is some mismodeling either in the IMP observa- 
tion model or in the GTDS data simulator. This is due to the large residual 
at the beginning of pass and to the presence of the monotonic decrease in the 
error with propagation only (i.e. , no filtering was done). There was a propa- 
gation of 20 minutes (spacecraft simulated time) prior to the first observation. 
However, a later test starting at the first observation showed little variation in 
the described results for the first few residuals. 

Further testing of the IMP ODS resolved the differences between the IMP model 
and the FORTRAN test model, but close scrutiny of the R&D GTDS simulated 
data model revealed an error there. The E&D GTDS simulated observation was 
modeled over an interval beginning 3/8 of a second later than the corresponding 
interval for the IMP model. The time displacement of the start in the R&D 
GTDS model was erroneous. This liappened when the correct first estimate of 
the observation was decreased by multiples of . 125 second to yield what is called 
the ambiguous observation. For the configuration of relay and target satellites 
considered here, three multiples of . 125 second are substracted from the range 
estimate internally. This ambiguous range is then used for calculating the 
modeling start time, thus accounting for the 3/8 second time error. This time 
discrepancy accounts for most of the difference (in either time or distance) 
between the final modeled ranges. The final difference is proportional to the 
current average range rate (i. e. , average time-rate of change of the range 
measurement). The correction of this error in the R&D GTDS model will be 
included in the next update. A corrected version is not available at this time for 
testing. 
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Runs wore nmdo in which differenco quotient upproxlmations to the range and 
range-rate partial derivatives with respect to the SMM (target satellite) state 
wore taken. The nunvorical approximations were made ns follows.* 

SR ehange in range 
SX^ change in X|,-component 


^^^^r . change in range rate (over the two observations after perturbations) 

SXj^ “ change in Xj^-convponent 

The results are shown in Table 2-5. Only a rough approximation of the partial 
dox’ivatlvos was expected to be obtained by these tests, because the difference 
quotients wore formed from quantities whoso changes were measured over 
10 -second Intervals. 

2.2 SYSTEM TESTS 

Three system tests of the application software wore performed. The first used 
IMP-16 simulated range sum and rango-sum-rato data as observations; the sec- 
ond used GT.DS simulated data; and tlie tliird used real preprocessed tracldng 
data. 

2,2.1 IMP-16 Simulated Data 

The IMP-16 ODS Kalman filter was unit tested by checking Its various compo- 
nent functions, as described in the previous sections. After the succesful 
completion of these tests, It was possible to run test cases of the entire orbit 
determination filter process. Those tost cases, which used observations gen- 
erated by the lMl^-16 ODS observation model, were designed to demonstrate 
that the filter is working. Tliis internal^ generated test data was preferred to 
tlie GTDS simulated tracking data because it presented no errors in the obser- 
vation model. 
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RANGE RESIDUALS (KILOMETERS) 
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Figure 2-5. Range Residuals (Propagation Without Filtering) 
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Table 2-5. Numerical Approximations of Observation 
Partial Derivatives for October 15 » 1979 


VARIATION 

NUMERICAU 

APPROXIMATION 

PARTIAL 
MODELED BY 
IMP 

RANGE (1 HOUR, 20 MINUTES, 0 SECONDS) 

■ 1 km 

-0.657 

-0.6559 

6y ■ 1 km 

0.750 

0-7507 

1 km 

-0.079 

-0.0789 

" 0.01 km/soc 

-0.300 

0.0 

» 0.01 km/soc 

0.100 

0.0 

6^ " 0.01 km/soc 

0.0 

0.0 

RANGE RATE (1 HOUR, 20 MINUTES, 10 SECONDS) 

5 “ 0,1 km 

-0.490 

-0.5658 

- 0.01 km/soc 

4893.05 

4920.18 
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Thd inothod of testing the filter NNlth internally' generated tost data is ns follows. 
The initial state of the target satellite (SMM) used in gonerat^' g the test data is 
perturbed relative to its conipononts by certain quantities. The filter is then 
started with tliis pertiu'bod target state and w^th the original starting state of the 
relay satellite (TORS) used for the data generation. The filter attempts to cor- 
rect the target state based on the internally' generated observations. The oidg- 
inal (unperturbed) orbit of the target satellite used during the data generation 
serves as the reference orbit for the target satellite (i.e. , for determining the 
position and velocity errors). 

All test cases whoso results are specified in this section used the same param- 
eters for filtering, tliese pammeters are shown in Table 2-6. The initial (iu\- 
porturbed) target satellite states are specified in Table 2-7. This test setup was 
repeated using both range and range rate (Run A), range only (Run B), and range 
rate only (Rim C) to correct the initially perturbed target element set. The 
results of these three rims are presented in Tables 2-8 and 2-9, 2-10 and 2-11, 
and 2-12 and 2-13, respectively. 

This set of initial target satellite state perturbations is particularly' difficult for 
the filter to handle. The first range observation i-esidual is very small (8 meters) 
compared to the rss error of the initial position perturbations (1732 meters). 
Thus, the input to the filter (i.e. , error signal) does not indicateoi large 
deviation from the expected state. Therefore, the state correction process is 
deferred until one or two observations later. It should be noted that for all 
test cases run with the given SMM and TDRS orbits, the tracking geometry 
(observability) is poor. This is confirmed by the very large (almost 1.0) 
correlations between X and Y and between X and Y. Such poor observability* 
makes tlie filtering process very difficult and unstable. 
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Table 2-6. Filter Test Parameters (Internally Generated Data) 


COMPONENT 

INITIAL STATE 
PERTURBATION 

INITIAL STATE 
COVARIANCE 

PROCESS NOISE 
COVARIANCE 

X 

1.0 km 

2 km^ 

1 e-6 km^ 

Y 

1.0 km 

2 km^ 

1 E-6 km^ 

Z 

1.0 km 

2 km^ 

1 E-6 km^ 

•k 

0,5 m/sac 

1 E-6 km^/sac^ 

1 E-12km^/sao^ 


0,5 m/sec 

1 E-6 km^/sec^ 

1 E-12km^/sec^ 

z 

0,5 m /sec 

1 E-6 km^/sec^ 

1 E-12 km^/sec^ 


NOTE: 

RANGE VARIANCE: 1 E-S KILOMETER SQUARED 
RANGE.RATE VARIANCE: 8 {CYCLES PER SECOND) SQUARED 
data rate {RANGE AND RANGE-RATE PAIR): ONCE EVERY 10 SECONDS UP TO 6 
MINUTES AND ONCE EVERY MINUTE AFTER 6 MiNUTES 
EDITING CRITERION: 48a 


Table 2-7. Initial (Unperturbed) Satellite States (at 1 Hour, 
20 Minutes, 00 Seconds on October 15, 1979) 


COMPONENT 

SMM 

TORS 

X 

2978.65815 km 

26306,1677 km 

Y 

-5942.30779 km 

-32582,1251 km 

Z 

1957.87311 km 

4797.49754 km 

X 

5.46081837 km/sec 

2.38132324 km/sec 

k 

3,90774318 km/sec 

1.94213659 km/sec 

z 

3.53571477 km/ssc 

0,133964630 km/sec 
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Table 2-8. Filter Test Results for Rim A (Range airi Range Rate) 


MINUTE 

POSITION ERROR Im) 

VELOCITY ERROR Im/sec) 

X 

mm 

2 

POSITION 

rss 

0 

X 

• 

Y 

t 

Z 

VELOCITY 

r« 

0 

1005 

994 

1001 

1732 

0.5 

0.5 

0.5 

0.866 

1 

1324 

1210 

658 

1878 

0.615 

0.270 

0.509 

0.843 


1326 

1194 

562 

1871 

0.528 

0.176 

0.471 

0,729 


1347 

1218 

804 

1986 

0.419 

0.073 

0.383 

0.572 


1295 

1137 

655 

1844 

0.371 

0.015 

0.374 

0.527 


922 

800 

272 

1261 

0.311 

0.127 

0.321 

0.465 


790 

667 

185 

1050 

0.308 

0.084 

0.322 

0,453 


534 

466 

42 

710 

0.255 

0,180 

0.251 

0.401 

8 

412 

348 

-11 

539 

0.232 

0.129 

0.204 

0.335 

9 

304 

274 

to 

1 

414 

0.106 

0.077 

0.087 

0,157 

10 


205 

-117 

329 

-0.004 

-0.011 

-0.020 

0.023 

11 


140 

-163 

280 

-0.090 

-0,101 

-0.005 

0.136 

12 


145 

-183 

284 

-0,185 

-0.145 

-0.032 

0.237 

13 


97 

-204 

264 

-0.210 

-0.189 

-0.014 

0,283 

14 

123 

126 

o 

CM 

1 

274 

-0.288 

-0.206 

-0.035 

0,356 

15 

101 

79 

-218 

253 

-0.286 

-0,225 

-0.007 

0,364 

16 

93 

118 

-199 

249 

-0.282 

-0.195 

-0.01 1 

0,343 

17 

80 

75 

-193 

222 

-0.249 

-0.192 

0,021 

0.315 

18 

61 

109 

-191 

228 

-0.289 

-0,198 

0.006 

0.350 

19 

57 

70 

-170 

192 

-0.227 

-0.174 

0.040 

0.289 

20 

53 

29 

-158 

169 

-0.197 

-0.164 

0.069 

0.265 


Table 2-9. State Error Covariance and Correlations for Rim A 



POSITION VARIANCE/CORRELATION (km^) 

VELOCITY VARIANCE/CORRELATION (km^/sec^i 

IVIKnU I c 

-1 


-I 

^XY 

4 

0^ 


PM 

5 

0.5446 

0,3856 

0.2060 

0.9989 

6,218 E *7 

4.525 E‘7 

9.367 E-7 

0,9413 

10 

0.1548 

8.285 E-2 

0.2049 

0.9989 

6.235 E-7 

3,449 E-7 

6.993 E-7 

0.9919 

15 

0.1595 

7.569 E-2 

0.3635 

0.9996 

2.328 E-7 

1,129 E-7 

4.010 E-7 

0.9895 

20 

0.1588 

6.876 E.2 

0.5335 

0.9997 

7.079 E-8 

3.777 E-8 

1.579 E-7 

0.9762 
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Table 2-10. Filter Test Results for Run B (Range Only) 


MINUTE 

POSITION ERROR |m) 

VELOCITY ERROR (m/jec) 

X 

Y 

Z 

POSITION 

rsi 

• 

X 

• 

Y 

• 

Z 

VELOCITY 

ru 

0 

1005 

994 

1001 

1732 

0,600 

0,500 

0.600 

0.866 

1 

1686 

1405 

236 

2132 

0.623 

0,301 

0.667 

0.888 


1419 

1272 

667 

1986 

0,620 

0.159 

0.484 

0.728 


1447 

1336 

922 

2176 

0,403 

0.016 

0.386 

0,658 


1336 

1170 

652 

1892 

0,355 

-0.032 

0,374 

0.517 


668 

666 

165 

958 

0.30B 

0.178 

0,293 

0,459 


696 

587 

124 

919 

0.300 

0.119 

0.300 

0.441 


519 

464 

33 

690 

0.234 

0.170 

0.233 

0.371 

8 

399 

336 

-31 

623 

0.178 

0.091 

0,161 

0.267 

9 

294 

264 

-107 

409 

-0,013 

-0.012 

-0.013 

0,022 

10 

227 

200 

-173 

349 

-0,161 

-0.129 

-0.057 

0.214 

11 

179 

137 

-223 

317 

-0,267 

-0.226 

-0.^)72 

0.360 

12 

159 

‘ 141 

-240 

321 

-0.339 

-0.267 

-0,088 

0,434 

13 

135 

94 

-244 

294 

-0.314 

-0.264 

-0.052 

0,414 

14 

119 

122 

-247 

300 

-0.380 

-0.271 

-0.067 

0,472 

15 

98 

76 

-242 

272 

-0,337 

-0,262 

-0,031 

0.428 

16 

90 

115 

-214 

269 

-0.312 

-0,215 

-0,031 

0,380 

17 

79 

74 

-197 

225 

-0,251 

-0,193 

0.006 

0.317 

18 

60 

108 

-195 

231 

-0,290 

-0.199 

-0.009 

0.352 

19 

55 

69 

-i6i:'i 

187 

-0,215 

-0,165 

0.029 

0,273 

20 * 

51 

28 

-153 

164 

-0,184 

-0.103 

0.067 

0.218 


Table 2-11. State Error Covariances and Correlations for Run B 



POSITION VARIANCE/CORRELATION (km^) 

i/ELOCITY VARIANCE/COHRELATION (km^/sec^! 

MINUTE 


a2 

2 

'’XY 

4 

4 

4 

'’X'f 

5 

0.4454 

0.3135 

0.1656 

0,9989 

6.190 S-7 

4.340 E-7 

9.338 E-7 

0,9522 

10 

0.1434 

7.608 E-2 

0.1950 

0,9991 

5.606 E-7 

3.080 E-7 

6.808 E-7 

0,9925 

15 

0.1577 

7.459 E-2 

0.3532 

0.9998 

1.463 E-7 

7.063 E-8 

3.967 E-7 

0,9851 

20 

0.1573 

6.791 E-2 

0.5297 

0.&'299 

4,191 E-8 

2.461 E-8 

1 ,546 E-7 

0.9675 
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Table 2-12, Filter Test Results for Run C (Range Rate Only) 


MINUTE 

POSITION ERROR (m) 

VELOCITY ERROR 

X 

Y 

Z 

POSITION 

rss 



• 

Z 

VELOCITY 

r» 


892 

2329 

891 


-0.730 

1,806 

0,296 

1.970 


2414 

1723 

2913 


-0,062 

0.274 

0.371 

0,465 


1227 

4658 

1237 


-1,214 

-2.822 

-0.240 

3.081 


1922 

1768 

2038 

3313 

-0,690 

-0.683 

0.156 

0.983 


3214 

1361 

3257 

4774 

-0.439 

-0.514 

-0,331 

0.763 


284 

1287 

834 

1560 

-0.796 

-1.177 

-0,228 

1.439 


1275 

1828 

1613 

2751 

-0,790 

-0.407 

-0.131 

0.898 

8 

1844 

2154 

2063 

3507 

-0.683 

-0.301 

-0.020 

0.747 

g 

1984 

2254 

2191 

3717 

-0,589 

-3.717 

-0.007 

0.928 

10 

1822 

2101 

2096 

3482 

-0,452 

-0.922 

0.007 

1.027 

11 

1516 

1805 

1901 

3028 

-0.316 

-0,985 

-0.020 

1.036 

12 

1192 

1517 

1684 

2561 

-0.268 

-0.952 

-0.116 

0.996 

13 

870 

1172 

1461 

2065 

-0.217 

-0.912 

-0.202 

0.959 

14 

596 

934 

1252 

1672 

-0.286 

-0.865 

-0.352 

0.977 

15 

349 

652 

1074 

1304 

-0.330 

-0.846 

-0.470 

1.023 

16 

136 

483 

901 

1031 

-0.484 

-0.844 

-0.639 

1.164 

17 

47 

270 

746 

794 

-0.587 

-0.869 

-0.769 

1.288 

18 

-214 

159 


652 

-0.798 

-0.895 

-0.918 

1.611 

19 

-345 

-1 


584 

-0.877 

-0.924 

-1.012 

1.627 

20 

-459 

-145 


600 

-0.980 

-0.952 

-1 .097 

1.752 


Table 2-13. State Error Covariances and Correlations for Run C 


MINUTE 

POSITION VARIANCE/CORRELATION (km^) 

VELOCITY VARIANCE/CORRELATION (km^/sec^l 

.,2 

“Y 

4 

'’XY 

a? 

“x 

4 



5 

1.541 

0.8403 

1.383 

0.4168 

6,430 E-7 

5.440 E-7 

1.179 E-6 

0,8589 

10 

0.9937 

0.5666 

0,9017 

0.1943 

7.719 E-7 

5.568 E-7 

1.709 E-6 

0.8646 

15 

’ 0.7933 

0.3599 

0.6228 

-8.909 E-2 

9.639 E-7 

5,078 E-7 

2.482 E-6 

0.9166 

20 

0,8475 

0.2678 

0.6342 

-0,1436 

1,195 E-6 

5.740 E-7 

3.352 E-6 

0.9602 
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2.2.2 GTDS Simulated Data 


Discrepancies between the IMP-16 observation model and the GTDS observation 
mo'Jel for SST data have been resolved. Test results will be supplied at a later 
date when the corrected model Is available in R&D GTDS. 

2,2.3 Preprocessed Tracking Data 

Real preprocessed tracking data will be available upon completion of the track- 
ing data preprocessor during Task 971. The results of filter tests using this 
data will be released when available. 
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SECTION 3 - SYSTEM TIMING AND PERFORMANCE 


One basic goal of the IMP-XC ODS development project was to determine the 
ability of the system to keep pace with Incoming tracking data in real time. 
Testing indicates that the computations for modeling and filtering one range and 
range-rate pair can be performed within 10 seconds (without drag), which is the 
expected data rate. A comprehensive evaluation must, of course, consider all 
input/output relevant to a pair of processed observations. 

3. 1 CURRENT SYSTEM TIMING WITH HEWLETT-PACKARD AND 
TEXAS INSTRUMENT TERMINALS 

The IMP- 16 processors used in the prototype orbit determination system dis- 
cussed here execute input/output under full central processing unit (CPU) con- 
trol only. Therefore, all input/output activity in a processor requires addi- 
tional lime for that processor. However, the design of the IMP- 16 ODS is such 
that the Data Base IMP (DBIMP) is responsible for input/output with the user, 
whereas the Computational IMP (COMPIMP) is responsible for the main compu- 
tational work of the orbit determination process (see the introductions to Chap- 
ters 5 and G of Reference 1). Tliis division of labor between the two px’ocessors 
is effective under continual processing of points. 

Table 3-1 indicates the approximate wall-clock times for the orbit propagator 
and filter reports for a range and range-rate pair under continual processing. 
The term "full output" applies to the present form of filter output, which uses 
one universal format for all output variable (array) types. The term "econo- 
mized output" refers to new formats to be implemented for filter reports. 

These new formats will save (printed) space and time. In addition, they will be 
more organized and, hence, more legible to the user. 

From Table 3-1 it may be seen that the EVIP-16 ODS can essentially keep pace 
with real time by using a fast (2400-band or faster) terminal. The Texas Instru- 
ment (TI) Silent 700, the terminal for the system, px’esents an upper physical 
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Table 3-1, System Speeds 


OUTPUT TYPE 

HP’ WALL-CLOCK TIMES (sec) 

T|2 WALL-CLOCK times (sec) 

FULL (PRESENT) 

12 

50 

ECONOMIZED (FUTURE) 

-10-12 

-35 


’HEWLETT-PACKARD 2640A INTELLIGENT TERMINAL AT 2400 BAUD 

^TEXAS INSTRUMENT SILENT 700 PRINTER TERMINAL AT 300 BAUD (HIGHER AVAILABLE BAUD 
RATES CANNOT BE USED DUE TO CONCENTRATION OF OUTPUT '.IMITATION OF PRINTING 
SPEED, AND LIMIT OF TERMINAL BUFFER MEMORY) 
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limit on system speeds of approximately 3,5 times real time while tracking, 
Depending upon the length of a data pass (nominally 20 minutes) and the rate at 
which observations are processed, the system should be able to catch up to 
i'eal time during data gaps (actual shadowing or Imposed gaps). For example, 
using a 20-second integrator step size and drag, the execution of a propagation 
thi'ough a 76-mlnute data gap (a 96-mlnute period minus 20 minutes of tracking 
data) takes approximately 11 minutes plus the Input/output time for periodic 
printout (plus 1-1/3 minutes for printouts produced once per minute). The 
approximate wall-clock time for filtering during the 20-minute pass would be 
approximately 3, 5 times 20, or 70 minutes. Thus, the total execution time for 
computation and input/output during the 96-mlnute period would be approximately 
88 minutes. 

3.2 TIMING OF INDIVIDUAL COMPONENTS 

The approximate breakdown of the computational time by component function is 
as follows; 

Time Required 


Component Function (seconds) 

Orbit propagator step without drag 3.3 

Orbit propagator step with drag 4.5 

Modeling range and range rate with partial^ 3.0 

Filtering range and range rate 4,0 


It should be noted that the times indicated in Table 3-1 may increase when the 
tracking data preprocessor or any other processor (in addition to the COMPIMP) 
is attached to the DPIMP for active communications. 


% 
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3.3 CURRENT/FUTURE CORE USAGE 

The brealcdown of core usage In the IMP CDS Is as follows j 



Core 

(bytes) 

Processor 

ramI 

PROM^ 

DBIMP 

5, OK 

6, OK 

COMPIMP 

C.OK 

15. OK 


With the addition of the tracking data preprocessor, core usage is expected to 
be IIK bytes of RAM and 9K bytes of PROM fcr DBIMP. COMPIMP core usage 
will remain the same. 

Current memory maps of DBIMP and COMPIMP are provided in Table 3-2, 


^Random Access Memory (not including base page). 
2 

Prograxnmable Read-Only Memory 
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Table 3-2 


Storage Maps 


COMPONENT 



MEMORY location’ 

RAM 

PROM 

DBIMP 

EXECUTIVE AND MESSAGE HANDLING SOFTWARE (DBMAIN) 

100-64C 

8000-8EAB 

UTILITIES (PTGTtI 

- 

8FOO-8F7C 

TERMINAL OUTPUT (TTYOUT) 

13BO-13F6 

9000-9469 

FLOATING-POINT PACKAGE (SFLPT) 


94F0-98AC 

OBSERVATION FILE (OBSFIL) 

660- 13 AO 

- 

COMPIMP 

COMPIMP EXECUTIVE AND MESSAGE HANDLER (COMPS) 

1280-142F 

8000-8220 

COMPIMP UTILITIES (UTIL) 

- 

8230-83FC 

FLOATING-POINT PACKAGE (FLPT) 

- 

8400-87CA 

MATH MODELS AND UTILITIES (DAGHST) 

180-25D 

8800-9170 

ORBIT PROPAGATOR (ORBIT) 

2B0-726 

9190-A7D0 

OBSERVATION MODEL (OBS) 

800-A24 

A800-B109 

SEQUENTIAL ESTIMATOR (KALFIL) 

21 00-23 A B 

B180-B918 

ORBIT DETERMINATION EXECUTIVE (COMP) 

1710-1979 

BA00-BC8E 


’addresses are hexadecimal. 
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